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Abstract: We study the propagation of ultra-high energy cosmic ray nuclei through 
the background of cosmic microwave and intergalactic infrared photons, using recent 
re-estimates for the density of the last ones. We perform a detailed Monte Carlo 
simulation to follow the disintegration histories of nuclei starting as Fe and reaching 
the Earth from extragalactic sources. We obtain the maximum energies of the arriving 
nuclear fragments as well as the mass composition as a function of the distance traveled. 
Cosmic rays with energies in excess of 2 x 10^° eV cannot originate from Fe nuclei 
produced in sources beyond 10 Mpc. 
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1. Introduction 

The cosmic ray (CR) spectrum is known to extend up to energies beyond 10^° eV, with 
the highest energy air showers observed having energies of 2-3x10^*^ eV. The origin 
and nature of these ultra-high energy (UHE) events are one of the pressing unsolved 
problems defying us today. The CR spectrum has the overall shape of a leg, and is well 
htted by power laws, whose index increases (spectrum steepening) for energies above 
the‘knee’ {E ~ 3x10^® eV), flattening again above the‘ankle’ (at E ~ 5x10^® eV). The 
CR composition becomes heavier for increasing energies around the knee, and the CR 
are probably mostly of galactic origin up to the ankle. Approaching the ankle, the CR 
composition seems to become lighter again, and the increasing rigidity of CRs does not 
allow anymore their conflnement into the Galaxy. Hence, CR fluxes are most probably 
of extragalactic origin above the ankle. There have been studies suggesting that the 
arrival direction of the highest energy events may be indicating that their origin lies 
in the local supercluster, but they are not conclusive. The small anisotropies observed 
may also be compatible with a cosmological origin of the highest energy events. 

The big difficulty which appears is that CR protons with E > 5 x 10^® eV, i.e. 
relativistic 7 factors > 5 x 10^°, are not able to propagate more than ~ 100 Mpc due 
to their energy losses by photopion production off the cosmic microwave background 
(CMB) photons, giving rise to the well known GZK cutoff |]^. At energies 2 x 10^° eV 
their mean free path is already only 30 Mpc. Heavy nuclei with smaller 7 factors, 
but comparable energies, also get attenuated but mainly by photodisintegrations off 
the intergalactic infrared (IR) background and off GMB photons, as well as by pair 
creation losses to a lesser extent Bill- 

A detailed study of the propagation of UHEGR Fe nuclei, including all the relevant 
energy loss mechanisms, was performed more than twenty years ago by Puget, Stecker 
and Bredekamp Q]. However, the estimates of the density of IR photons employed 
then were about an order of magnitude larger than the new empirically based estimates 
obtained using the measured emissivity of IRAS galaxies |^. In the light of the lower 




IR background densities inferred recently, it was suggested that UHECR nuclei may 
propagate much longer distances unattenuated [0], so that the events with energies 
2-3x10^° eV could have possibly originated as Fe nuclei produced at distances up 
to 100 Mpc, and hence in particular in the whole local supercluster^. However, as we 
showed in a recent letter p, at energies larger than 10^° eV it is photodisintegration off 
CMB photons (and not off IR ones) which dominates the opacity for Fe disintegration. 
This implies that the maximum energies with which the surviving fragments can reach 
the Farth are not signihcantly changed (for distances below a few hundred Mpc) with 
the new estimates of the IR density. In particular, for sources at distances of 100 Mpc 
the maximum energies of the surviving fragments do not exceed ~ 10^° eV, and can 
arrive to 2 x 10^° eV only for distances below 10 Mpc. 

The aim of the present paper is to re-evaluate in detail the propagation of heavy 
nuclei, following the photodisintegration histories by means of a Monte Carlo which 
includes all relevant processes, much in the spirit of the original Puget et al. paper [^]. 
From this we can establish all the effects resulting from the new estimates of the IR 
background density. In particular, we obtain the hnal mass composition and energy as 
a function of the distance to the source, as well as the possible fluctuations in these 
quantities which may arise from the particular way in which the photodisintegration 
takes place in each case. We also study the effects of pair creation losses, which turn 
out to be important in some cases for the determination of the hnal mass composition. 


2. The propagation of heavy nuclei 

As we said before, CR with energies above the ankle are most probably extragalactic. 
This means that in their journey they may be attenuated by the interactions with the 
photon background. This background consists essentially of the microwave photons 
of the 2.7°K cosmic background radiation and, at larger energies, of the intergalactic 
background of IR photons emitted by galaxies. The background of optical radiation 
turns out to be of no relevance for UHECR propagation. 

Although the CMB density is well known, the intergalactic IR one cannot be mea¬ 
sured directly and has to be estimated from the observation of the spatial distribution, 
IR spectra and emissivity of the galaxies which are sources for this IR emission. This 
was done recently by Malkan and Stecker 0, who obtained a result which is about 
an order of magnitude smaller than previous estimates. We will then adopt for the 
spectral density of the IR background 

^=l.lxl0-hAj"“cm-’eV-‘ (2.1) 

de VeV/ 

^It has to be stressed that Fe nuclei are good candidates for UHECRs, due to their high abundance 
in supernova environments and their large value of Z^ which enhances the energy achievable in the 
acceleration process. 



for photon energies e in the range between 2 x 10“^ eV and 0.8 eV. This is a factor of 10 
smaller than the “high infrared (HIR)” density adopted in ref. [Q, and is in the upper 
range of the recent estimates. In order to quantify the possible effects of an optical 
intergalactic background, we just modeled this last with a Planckian distribution with 
T = 5000°K and a dilution factor of 1.2 x 10“^®, as in |Q. 

UHECR nuclei propagating through these photon backgrounds will loose energy 
mainly by two processes: 

i) photopair production, which has a threshold corresponding to photon energies in 
the rest frame of the nucleus of 2meC^ ~ 1 MeV. This process was studied in detail by 
Blumenthal 0 , and the main contribution arises from interactions with CMB photons. 


We used for the energy loss rate the expressions given in ref. 


ii) photodisintegration losses, for which the rate of emission of i nucleons from a 
nucleus of mass A (with cross section j) is given by 




1 

2 ’ 7 ^ 


'0 


- de dn /■27^ 


( 2 . 2 ) 


where 7 is the relativistic factor of the nucleus, e the photon energy in the observer’s 
system and e' its energy in the rest frame of the CR nucleus. 

The cross sections for photodisintegration contain essentially two regimes. 

At e' < 30 MeV there is the domain of the giant resonance and the disintegration 
proceeds mainly by the emission of one or two nucleons. At higher energies, the cross 
section is dominated by multi-nucleon emission for heavy nuclei and is approximately 
flat up to e' ~ 150 MeV. We fitted the various aA,i with the parameters in Table I and 
II of ref. [Q. A useful quantity to estimate the energy loss rate by photodisintegration 
is given by the effective rate 

A A 

Reff,A = (2-3) 


Since neglecting pair creation processes one has that photodisintegrations alone 
lead to E~^dE/dt = A~^dA/dt, the energy loss time for photodisintegration is then 
A/Reff^A- The different contributions to this quantity are plotted in Fig. 1. We show 
separately the contributions to the disintegration from CMB, IR and optical photons 
for Fe nuclei, together with the total one (solid line) and the photopair creation energy 
loss rate^. 

It is apparent that the optical background has no relevant effect, that the IR one 
dominates the photodisintegration processes below 10^° eV and the CMB dominates 
above 10^° eV. The pair creation rate is relevant for Fe energies 4 x 10^® eV-2 x 10^° eV 
(i.e. 7 factors ~ 1-4 x 10®), for which the typical CMB photon energy in the rest frame 
of the nucleus is above threshold (> 1 MeV) but still well below the peak of the giant 
resonance (~ 10-20 MeV). The effect of pair creation losses is to reduce the 7 factor 
of the nucleus, obviously leaving A unchanged. 


^Looking at this and the following figures, it is important to keep in mind that 1 Mpc= 1.03 x 10^"* s. 







Figure 1: Effective energy loss time for Fe photodisintegration off microwave (CMB), infrared 
(IR) and optical (OP) photons, as well as the total one (solid line) and the pair creation (PC) 
energy loss time. 

3. Results 

Using the rates just discussed, we performed a Monte Carlo simulation in order to follow 
the possible disintegration histories of Fe nuclei. In hgures 2 and 3 we plot the hnal 
mass (i.e. A) and energy E of the heaviest fragment surviving from the disintegration 
process. We show in these hgures the results of simulations with initial values of the 
relativistic factor 70 = 1 x 10^°, 4 x 10® and 1 x 10®. The curves shown for each value of 
7 o correspond, in Fig. 2, to the average value (A) from all the simulations (solid line) 
and the region (between the two dashed lines) including 95% of the simulations^. This 
gives a clear idea of the range of values which can result from huctuations from the 
average behaviour. 

To further understand the relevance of the different processes and the impact of 
the new determinations of the IR density, we also plot the results for {%) obtained in 
simulations which do not include pair creation processes (dotted lines) and also the 
results we would obtain (dot-dashed line) with an IR density a factor of ten larger (i.e. 
the HIR density of ref. [§]). Figure 3 is similar but for the hnal values of the energy. 

From these hgures we can draw the following conclusions: 

i) For large initial energies {Eq > 2.5 x 10^® eV, i.e. 70 > 5 x 10®), both the ehects 
of the IR photons and of pair creation processes are of no relevance along the whole 
journey of the nucleus, and the energy losses are essentially due to photodisintegration 
oh CMB photons alone. 

ii) At 7 o < 5 X 10® the pair creation losses start to be relevant, reducing the value of 
7 signihcantly as the nucleus propagates distances 0(100 Mpc). The ehect is maximum 


^Only 2.5% of the simulations are below the lower curves and 2.5% are above the upper ones. 
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Figure 2: Evolution of the mass number A of the heaviest fragment surviving photodisinte¬ 
gration vs. travel time t. The initial 7 factors considered are 1 x 10^® (I.EIO curve), 4 x 10® 
(4.E9 curve) and 1 x 10® (1.E9). 
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Figure 3: Similar as Figure 2 but for the final energy E vs. t. 

for 7 o ~ 4 X 10® but becomes small again for 70 < 1 x 10®, for which appreciable effects 
only appear for cosmological distances (> 10^ Mpc), as can be simply understood from 
Fig. 1. The effect of neglecting pair creation losses translates into keeping 7 = 70 
constant during the propagation, and this enhances the photodisintegration rates and 
then reduces (A) more rapidly. 

in) Also for 70 < 5 x 10® the reduction in the IR density adopted has sizeable 
effects. In this respect point ii) is relevant, since pair creation losses shift the valnes of 
7 towards a domain where IR photons become increasingly important with respect to 
CMB ones. With the new valnes of the IR density the effects of photodisintegrations 
























become small already for 70 ~ 1 x 10 ® if we consider propagation distances below 
10^ Mpc (i.e. for t < 10^^ s). 

iv) The effects of neglecting pair creation losses are less pro nonneed in Fig. 3. For 
instance, for 70 = 4 x 10® the average energies with and withont pair creation processes 
are similar np to t ~ 10^® s while the (A) valnes differ sizeably already for t ~ 3 x 10^® s. 
This is dne to a partial cancellation between the effects of the evolntion of 7 and of A 
in the valnes of the hnal energy {E = rUp'yA), since neglecting pair creation losses does 
not allow 7 to decrease bnt makes instead A to drop faster^. 

v) The effects of flnctnations dne to different photodisintegration histories are not 
negligible. They give a spread in A (and E) of the order of 10% (considering the 95% 
probability range) for (%) ~ 40 bnt relatively larger for smaller (%), since variations 
AA ~ 10-15 at a given time t can appear between different simnlations. 



Figure 4: (A) vs. t for more sample values of 70 . 


In hgnres 4 and 5 we plot more sample valnes of {%) and {E), for valnes of 70 = 
2 X 10^®, 1 X 10^® and ( 8 , 6 , 4, 2, l)xl0®. Looking at Fig. 5 it is easy to infer the 
maximnm energies which can be obtained from Fe nnclei injected at any hxed distance 
d. In particnlar, for d = 100 Mpc {t = 10^® s) the maximnm average energy is E^ax — 
8 X 10^® eV and originates from 70 ~ 2-4x 10®. Comparing with Fig. 4 we see that these 
maximnm energy events wonld correspond to fragments with masses A{Emax) — 30-50, 
i.e. a rather heavy composition. 

Flnctnations from the average behavionr can only slightly increase the maximnm 
attainable energies, and this is illustrated with the dashed line, which represents the 
upper boundary of the 95% CL ranges (i.e. 97.5% of the simulations are below this 
curve) for any initial value of 70 . 

'^This in particular shows that the inclusion of pair creation losses does not modify the maximum 
attainable energies computed in |^. 






For source distances d = 10 Mpc, average energies up to E^ax — 2 x 10^° eV can 
result, and an interesting pile-up effect is observed since a broad range of initial energies 
(with 7 o ~ 4-8 X 10®) lead to approximately the same hnal energy ( Emax)- This can 
produce a bump in the spectrum from sources at these distances if indeed the highest 
energy events originate from heavy nuclei. Due to the spread in values of 70 at E^ax, 
we see from a comparison with Fig. 4 that there will also be a wide spread in the hnal 
composition, with A{Emax) — 10-45. Events with energies 2-3 x 10^® eV may appear 
as low probability huctuations from the mean behaviour if d ~ 5-8, having initially 
7 o > 10^® and a low mass hnal composition [A < 10). 
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Figure 5: {E) vs. t for several sample values of 70 . The dotted line indicates the upper 
boundaries of the 95% probability range obtained in the simulations. 

For smaller source distances (d ~ few Mpc), events with E ~ 3 x 10^® eV could be 
heavier nuclei with smaller initial energies (70 > 6 x 10 ®). 

For very large values of 70 (70 > 2 x 10^®), the heavy nuclei completely disintegrate 
in less than 10 Mpc, and the photopion production (not included here) becomes the 
main attenuation process for the secondary nucleons, which are then subject to the 
usual GZK cutoh. 

In conclusion, the main implication of the lower values of the IR density recently 
estimated is to increase the mean free path of the heavy nuclei with initial 7 factors 
below ~ 5 X 10®, for which most CMB photons are below the peak of the giant resonance 
for photodisintegration. Due to the fragmentation of the nuclei by photodisintegration 
and the pair creation energy losses, the hnal energies of the fragments are typically 
below 2 X 10^® eV for travel distances ~ 10 Mpc, and below 10^® eV for distances 
~ 100 Mpc. The new value of the IR density is then of little help in the attempts to 
understand the highest energy events observed {E ~ 2-3 x 10^® eV), which could not 
have originated as heavy nuclei at distances beyond ~ 10 Mpc. The lack of obvious 




candidate sources at closer distances p leave the nature and origin of these events still 
a mistery. 
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